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Abstract. We consider the propagation of wave packets for a one-dimensional 
nonlinear Schrodinger equation with a matrix-valued potential, in the semi- 
classical limit. For an initial coherent state polarized along some eigenvector, 
we prove that the nonlinear evolution preserves the separation of modes, in a 
scaling such that nonlinear effects are critical (the envelope equation is non- 
linear). The proof relies on a fine geometric analysis of the role of spectral 
projectors, which is compatible with the treatment of nonlinearities. We also 
prove a nonlinear superposition principle for these adiabatic wave packets. 



We consider the semi-classical limit e — >■ for the nonlinear Schrodinger equation 



where A g R. The data ipQ and the solution tp'^{t) are vectors of C^, N ^ 1. The 
quantity |'!/''^lc« denotes the square of the Hcrmitian norm in of the vector 
Finally, the potential V is smooth and valued in the set of by Hermitian 
matrices. Such systems appear in the modelling of Bose-Einstein condensate (see [T] 
and references therein). 

Definition 1.1. We say that a function f is at most quadratic if f E C°°(R) and 
for all k^2, /W e L°°(R). 

We make the following assumptions on the potential V: 

Assumption 1.2. (1) We haveV{x) = D{x) + W{x) with D,W e C°° {11,^1^ '"^), 
D diagonal with at most quadratic coefficients, and W symmetric and bounded as 
well as its derivatives, W G VF°°'°°(R). 

(2) The matrix V has P distinct, at most quadratic, eigenvalues Xi, . . . , Xp and 
(1.2) 3co,no e R+, Vj ^ fc, Vx e R, |A,(x) - Afc(x)| ;? co (a;)""" . 

Under these assumptions (the first point suffices) , we can prove global existence 
of the solution ip^ for fixed e > 0: 

Lemma 1.3. IfV satisfies A ssumvtion \1.2\ and ipo G ^^(R); there exists a unique, 
global, solution to (jl.ip 
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The L^-norm of ip'^ does not depend on time: ||'i/''^(^)l|L2(R) = ||V'olU^(R)i ^ 
The proof of this lemma is sketched in Appendix [A] 

In this nonhnear setting, the size of the initial data is crucial. As in [4], we 
choose to consider initial data of order 1 (in L^), and to introduce a dependence 
upon s in the coupling constant (note that the nonlinearity is homogeneous) . This 
leads to the equation 

and we choose the exponent /? = 3/4, which is critical for the type of initial data 
we want to consider (coherent state) when the potential V is scalar (see ;4j). We 
are left with the nonlinear semi-classical Schrodinger equation 

(1.3) ^e^tr + ^^lr = Vix)r+Ae'/'\r\h«^' ; ^\t=o = ^o- 
We focus on initial data which are perturbation of wave packets 

(1.4) V§(^) = e-i/4e^«"(-~-«)/^a (^^) x(^) + r'oi^), 
where the initial error satisfies 

(1-5) IkolU^CR) + \\xro\\L2(R) + \\edxro\\L2{-R) = 0{e'^) for some k > ^• 

The profile a belongs to the Schwartz class, a e iS(R; C), and the initial datum is 
polarized along an eigenvector xi^) G C°°(R; C^): 

V{x)x{x) = Xi{x)x{x), with \x{x)\c« = 1- 
Note that Ai is simply a notation for some eigenvalue, up to a renumbering of 
eigenvalues. The L^-norm of ipQ is independent of e, ||V'oIIl2(r) = ||a||L2(R). As 
pointed out above, this is equivalent to considering (jl.ip with initial data of the 
same form (|1.4p . but of order e^/"' in L^(R). The evolution of such data when a is 
a Gaussian has been extensively studied by G. Hagedorn on the one hand, and by 
G. Hagedorn and A. Joye on the other hand, in the linear context A = (see [HIS])- 
These data are also particularly interesting for numerics (see [l^ and the references 
therein). 

Because of the gap condition, the matrix V has smooth eigenvalues and eigenpro- 
jectors (see [11 J. Besides, the gap condition (|1.2p also implies that we control the 
growth of the eigenprojectors (see Lemma lC.2p . Note however than in dimension 1 
(x G R), one can have smooth eigenprojectors without any gap condition. We 
explain this fact below and give an example of projectors that we can consider; we 
also illustrate why things may be more complicated in higher dimensions {d ^ 2). 

Example 1.4. For N — 2 and x e R, consider 

(1.6) Vix) = {ax^ + b)ld + ( "5"^} ''^';\ ) , 

for a,b £ R, and u and v smooth and bounded with bounded derivatives. Such a 
potential satisfies Assumption 11.21 Its eigenvalues are the two functions 

{x) =ax'^ + b± ^yu{x)'^ + v{x)'^. 

These functions are clearly smooth outside the set of points xq such that w(xo)^ + 
v{xo)'^ — 0. Besides, for such points, one can renumber the modes in order to 
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build smooth eigenvalues. More precisely, observe first that if u{x)'^ + — 
0{{x — xq)^) close to xo, the functions are smooth close to xq. Moreover, if 
u{x)'^ + v{x)'^ = {x — xqY J{x) with /(a;o) 7^ 0, necessarily /(xq) > and k — 2p, 
so we have 

X^{x) = ax^ + b±\x- xo\P^/f(x). 

For p even these functions are again smooth. However, when p is odd, they are no 
longer smooth and we perform a renumbering of the eigenfunctions, observing that 

X I— >■ ax^ + b + {x ~ xqY \/ f {x) 

are smooth eigenvalues of V close to xq . 

Example 1.5. Resume the above example, with now x e R'', d ^ 2. The smoothness 
of the eigenvalues is no longer guaranteed: suppose u{x) = xi and v{x) = X2, then 
the functions A± are not smooth and one cannot find any renumbering which makes 
them smooth. 

Example 1.6. For an example of a potential which satisfies II. 2[ we simply choose 
V as in (11.61) with 



(1.10) idtu+-dlu^-\'i{x{t))y''u + k\u\^u ■ u{0,y)^a{y), 



Cuuix) = Cyv{x) = (x) "° , cl + cl ^ 0. 
1.1. The ansatz. We consider the classical trajectories {x{t),£^{t)) solutions to 

(1.7) x{t) = at), iit) = -VAi(x(t)), xiO) = XQ, m - Co. 

Because Ai is at most quadratic, the classical trajectories grow at most exponen- 
tially in time (see e.g. [4]): 

(1.8) 30 0, |e(i)l + Kt)l<e^*. 
We denote by 5* the action associated with {x{t),(^{t)) 

(1.9) ^(i)=£(^le(-s)p-Ai(x(.))) ds. 
We consider the function u — u{t,y) solution to 

and we denote by (p^ the function associated with u, x, ^, S by: 

(1.11) ^'{t,x) = S-'/\ (^t, ^^1^) 

Global existence of u and control of its derivatives and momenta are proved in [5] . 
More precisely, we have the following result. 

Theorem 1.7 (From [3]). Suppose a G iS(R). There exists a unique, global so- 
lution u e C(R;L2(r)) n Lf„^(R;i'^(R)) to (frTO| . In addition, for all k,p e N, 
{y)'" dPueC{R;L^{Il)) and 

(1.12) Vfc,peN, 30Q, VteR+, IKy)' ^^^/(t, OIIl^ir) < e^*. 

In particular, note that dlyu{t,-) is in L°° for all p e N. These results have 
consequences on Lp^ . As far as the L°° norm is concerned, we infer, using 

(1.13) VpeN, ||(£a,)V(0||Loo <£-i/W. 
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We use the time-dependent eigenvectors constructed in [8 (see also [9, and ,14i). 
To make the notations precise, we denote by dj the multipHcity of the eigenvalue Xj, 
1 =^ j < P (note that Ei^j^jP '^J = 

Proposition 1.8. There exists a smooth orthonormal family (x^(^: 2;)) ^^^^^ such 
that for allt, (x^(t, a;)) ^^^^^^^ spans the eigenspace associated to Ai, x^i^:^) — x{^) 
and for m € {1, • • • , di}, 

(1.14) {x"'{t,x),dtx'{t,x) + md.x'{t,^))c« = ^- 

Moreover, for £ e {1, . . . , di}, k,p G N, there exists a constant C — C{p, k) such 

that 

\dfd'y{t,x)\^,^Ce^' (^)('=+p)(i+"o)^ 
where ng appears in (jl.2l) . 

Note that equation (|1.14p for m = is true as soon as the eigenvector x^ is 
normalized and real-valued. 

Equation (|1.14p is often referred to as parallel transport. These time-dependent 
eigenvectors are commonly used in adiabatic theory and are connected with the 
Berry phase (see [H]). Their construction is recalled in Section[2l where the control 
of their growth is also established. 

Notation. In the case of a single coherent state, we complete the family (x^{t, x)) 
as an orthonormal basis (x^) i^jssp of as follows: 

e e 

vjA^Kdj span; 

1.2. The results. We prove that there is adiabatic decoupling for the solution 
of (|1.3p with initial data which are coherent states of the form (|1.4p : the solution 
keeps the same form and remains in the same eigenspace. 

Theorem 1.9. Le< a G 5(R) and rf, satisfying p.5p . Under Assumption ll.Sl con- 
sider ■0'^ solution to the Cauchy problem (|1.3p - (|1.4p . and the approximate solution 
iff^ given by (jl.lip . There exists a constant C > such that the function 

w%t,x) = 'ip%t,x) - (p''{t,x)x^{t,x), 

where x^ is given by Provosition \1.8l satisfies 

sup i\\w%t)\\L2 + \\xw%t)\\L2 + \\sd.,w'it)\\L2)^0. 
ItKCloglogi 



• For j ^ 2 and 1 ^ £ ^ dj, Xj — Xji^) does not depend on time, 

• For j ^ 2, {Xj)i^j^dj spans the eigenspace associated to Xj. 



This adiabatic decoupling between the modes is well-known in the linear setting 
and is at the basis of numerous results on semi-classical Schrodinger operator with 
matrix-valued potential in the framework of Born-Oppenheimer approximation for 
molecular dynamics. On this subject, the reader can consult the article of H. Spohn 
and S. Teufel [13] or the book of S. Teufel |T3] for a review on the topic (see also [1] 
for an adiabatic result in a nonlinear context and |10) for application of adiabatic 
theory to the obtention of resolvent estimates). 

Remark 1.10. Suppose that V depends on e with = D + eW, where D and W 
are as in Assumption II. 2 1 this is so in the model presented in [1]. Then the above 
result remains true for \t\ ^ Clog (-): we gain one logarithm. See Remark |4. II for 
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the key arguments. Also, the assumption on the initial error can be relaxed: to 
prove the analogue of Theorem 11.91 with an approximation in up to Clogl/e, 
(|1.5p can be replaced with 



In contrast with the general framework of this paper, no rate is needed: the rate in 
(jl.Sp is due to the fact that we cannot use Strichartz estimates here. 

It is also interesting to analyze the evolution of solution associated with data 
which are the superposition of two data of the studied form. We suppose 



where both functions tpi and i/?! have the form (|1.11[) . for two eigenvectors of V, 
Xi and X2, and phase space points (a;i,Ci) and (a;2,^2)- We assume 

{Xi,xi,£,i) ^ (X2,a;2,6) ■ 
We associate with the phase space points {xj,^j), j E {1,2} the classical trajectories 
{xj{t),(^j{t)), and the action Sj{t) associated with Xj such that 

V{x)xjix) = \jix)xjix). 

Note that we may have Ai = A2. Let us denote by x1(0 i^j^P a time-dependent 

orthonormal basis of eigenvectors defined according to Proposition 12.11 (see also 
Proposition II .81 above) with Xi(0,a;) = Xi{x), X2ix) = Xi(0,a;) if Ai — A2, X2{x) = 
X2(0,a;) otherwise, and by </?^ the ansatz defined by (jl.lip . To unify the presenta- 
tion, we write 



Note that if Ai ~ A2, one recovers the condition Ei 7^ E2 of [4 . The proof of 
Theorem 11.111 follows the same lines as in [U Section 6]. The constant F controls 
the frequencies of time interval where trajectories cross. 

Remark 1.12. In finite time, the situation is different whether Ai = A2 or not. If 
Ai = A2, the superposition holds in finite time without any condition on F; this 
comes from the fact that the trajectories xi{t) and X2{t) only cross on isolated 
points (see (5). However, if Ai ^ A2 one may have xi{t) = X2{t) on intervals of 
non-empty interior: the condition F ^ prevents this situation from happening. 
For example, if 



lkolU=(R) ^ as e ^ 0. 



i'oix) = V'l (0, x)xi{x) + (/3|(0, x)x2{x), 




Theorem 1.11. Set Ej = ^- + Xj{xj) for j e {1, 2} and suppose 
F=inf Xi{x) - X2{x) ~ {El - E2) >0. 




0. 
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with V smooth and at most quadratic, we have \i{x) = — 1 and A2(a;) — v[x) + l: 
classical trajectories for both modes, issued from the same point of the phase space, 
are equal. 

1.3. Strategy of the proof of Theorem 11.91 The proof is more complicated 
than in the scalar case j4j, due to the fact that the spectral projectors do not 
commute with the Laplace operator. From this perspective, a much finer geometric 
understanding is needed and we revisit [H [71 HI [131 [S] by adapting to our nonlinear 
context ideas contained therein. 

Observe first that the function Lp^ satisfies 

(1.15) ledt^' + ^92^^ ^ %{t,x)v' +Ke^'^\v%.v', 
where 

Te{t,x) = Xi{x{t)) + K{x{t)){x - x{t)) + ]^\"{x{t)){x - x{t))\ 

This term corresponds to the beginning of the Taylor expansion of Ai about x{t). 
Therefore, the function a;) — ip^{t, x) — x)x^it, x) satisfies w[(^q = 

iedtw,{t,x) + —dlw''{t,x) - V{x)w^{t,x) = eNL {t, x) + eL^t, x) 

where 

= zdtx'^' + ed.x'd.ip' + ^^'dlx' + (Ai(x) - %) ^'x'- 
Since tp'^ is concentrated near x — x{t) at scale ^/e, we have 

(Ai(x) - %)ip' - O (£3/2eC^*) in L^K), 

where we have used Theorem 11.71 The term a priori presents an 0(1) con- 
tribution, which is an obstruction to infer that is small by applying Gronwall 
Lemma. Observing that in view of the estimates on the classical fiow (see (jl.Sp ) 

ed^if" = i^iijif" + O (v^e") in ^^(R), 

we write, 

= ^ {dtx' + md.x') V' + O (V^c^') in L^(R). 
The choice of the time-dependent eigenvectors ensures that for all time, the 0(1) 
contribution of U is orthogonal to the first mode (the eigenspace associated with 
A^). Then, to get rid of these terms, we introduce a correction term to w'^ . We set 

e'{t,x)^w'{t,x)+eg'{t,x), g'it,x) = ^ ^ gl^it,x)x%x), 

where for j ^ 2 and for 1 ^ £ ^ dj, the function gj^{t,x) solves the scalar 
Schrodinger equation 

£2 

(1.16) iedtgle + y^^^^ ~ >^ji^)9j,i = ; 9j,i\t=o = 0' 
where 

(1.17) r,4t,x) = -I {dtx\t,x) + i{t)dxx\t.x) , X,'(a;))^„ . 
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The function 9'^{t) then solves 



(1.18) 



with 



iedt6'{t,x) + —die" it, x) = Vix)0'{t,x) + eNL''{t,x) + eL'{t,x), 



t=0 — 'Oi 



(1.19) NL- = A£i/2 (i^e^i +0e^ sg%.{^'x' +0'- ef) - l^lV^x') , 
(1.20) 



= L'+[ iedt + —dl - Y{x) ) g^{t, x) 



r ^2 



where the ©(i/ee*-^*) holds in L^. The proof of the theorem then follows from a 
precise control of the functions x'j ^nd g^^, which is achieved in Sections [2] and [3l 
respectively. Then, the analysis of 9^ as e goes to zero by an energy method is 
presented in Section ID 



2. The family of time-dependent eigenvectors 

In this section we prove Proposition ll.Si recalling the construction of the eigen- 
vectors satisfying (|1.14p . and analyzing the behavior of their derivatives for large 
time. We follow the proof of [8]. More generally, we prove the following result which 
implies Proposition ll.Si We consider the Hamiltonian curves of i|CP + Aj(a;), that 
we denote by {xj{t),£_j{t)). 

Proposition 2.1. There exists a smooth orthonormal basis of (x^(t,x)) i^e^dj 

such that for all t, {Xj{t,x)^ spans the eigenspace associated to Xj, with 

X^{0, x) = x{^) for m e {1, • • • , dj}, 

(X7(t, x),dtx% x) + ij{t)d^X% a;))^„ = 0. 

Moreover, for £ G {1, . . . ,dj}, k,p G N, there exists a constant C = C{p,k) such 
that 

|9ra,^x^(i,x)|^„^Ce^*(x)('=+^)(i+"«), 

where tiq appears in (11.21) . 

Proof of Proposition \2.1\ We consider a smooth basis of eigenvectors (x^ (0)) i^i^dj 

such that Xi(0) = X- Then, we denote by T{j{x) the smooth eigenprojector associ- 
ated with the eigenvalue Xj {x) and define 

Kj{x) ^ ~i[W,{x),d^Wj{x)] . 

We set z — x — Xj (t) and we consider the Schrodinger type equation 

(2.1) idtY^it, z) - £,j{t)Kj{z + x,{t))Yf{t, z) ; r/(0, z) = x'ji^M + z). 
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Let US prove that the vector Yj{t, z) is m the eigenspace of Xj{xj{t) + z). Indeed, 
the evolution of z) = (Id - Tlj{xj{t) + z)) Yf{t, z) obeys to Zj(0, z) = and 

- W(id - n,(x,(t) + z))[n,{x,{t) + z), d,n,{xj{t) + z)]Yf 
= -^,{t)d.n,{x,{t) + z){id - n,{xj{t) + z))Y^' 

= -^,(t)a,n,(a;j(t) + z)Zj 
where we have used dx^j — dx(J^) = lijdx^j + {dx^j) Hj, whence 

n,(a,n,)n, = n, (n,a,n, + (a,n,)n,)n, = 2n,(a,n,)n, = o. 

Therefore, Zj{t) satisfies an equation of the form 9tZj = z)Z^^, which combined 
with ^|(0) = 0, imphes Z^^it) — for all f € R: the vectors Yj(t, z) are eigenvectors 
of V{xj{t) + z) for the eigenvalue Xj{xj{t) + z). 

Besides, since ^j{t)Kj{z + Xj{t)) is self-adjoint, Y-'{t, z) is normalized for all i, and 
the family {Y^)i<^(^<^d. is orthonormal. We define by 

(2.2) x%x)^Y^{t,x^x,{t)) 

and we obtain an orthonormal basis of eigenvectors of V{x). 
It remains to check that (I1.14p holds. We have 

dtX^+^dxx] - dtY^'{t,x-x,{t)) 

= zi,{t)K,{x)x] 

- -i,{m,{x),dxli,{x)]x], 

whence 



= -e,(i)(n,[n,,9,n,]n,x,U.r 



since xj'"^ = ^jX'j''^- We then observe that 11^ = Ilj implies 

n,[nj,a,nj]nj - n^^^n^n^ - n.a.n.n^ = o. 

This concludes the first part of Proposition 11.81 It remains to study the behavior 
at infinity of the vectors Xj(^i 2;) and of their derivatives. 

By the definition of in (|2.2I) . it is enough to prove 

\dfd^xY,'it,z)\c^ < e^*(x,(t) + z)(^+'=)(i+"«). 

For this, we crucially use the estimates of Lemma IC.21 and we argue by induction. 
Let us first consider the case p = 1 and k — 0. By Lemma FC. 21 we have |iirj(x)| < 
{x)^~^"° , whence (|2.1I) gives 

\d,Ylit, X - x,m < \^,mKjix)\ < e^* (a;)^-*-"" . 

Let us now suppose k ^ 1 and p = 0. We observe that d^Yj{t, z) solves 

(2.3) 



idtd'^YHt, z) = ^ii,{t)K,{z + x,{t))dlYHt, z) + /(i, z), 



^0 V 

9^^r,(0,z) = 9^X,'(0,2 + a;,(0)), 



A NONLINEAR ADIABATIC THEOREM FOR, COHERENT STATES 



9 



where 

f{t,z)= c^^,{t)d]K,iz + x,{t))d^-'Yj{t,z), 

for some complex numbers independent of t and z. We obtain 

d^YHt,z)=U,{t,0)d^^X-{z + Xj{0))+ fu,{t,s)f{s,z)ds, 



where t{j{t, s) denotes the unitary propagator associated to p.ip (when the initial 
time is equal to s). We have by Lemma rC.2l 



\d^Y^'{0,z)\c^ = |5,^x^(0,z + a:,(0))|^„ < {z + x,{0))'^'+-'^\ 

therefore the induction assumption 

V7 e {0, . . . , fc - 1}, \d2Ylit, z)|c« < e" {xj{t) + 

implies, along with Lemma IC. 21 

\d'.Yf{t,z)\a.<e^Hx,it) + z)'^'^-«\ 

We have obtained the estimate for p — 0, k E N, and for p ~ 1, k — 0. Note that 
Equation ()2.3|1 yields 

and allows to prove the general estimate for time derivatives by an induction ar- 
gument which crucially uses the fact that we have an exponential control of the 
derivatives in time of £,j{t). This property follows by induction from (|1.7I) . (|1.8I) . 
and the fact that Xj is at most quadratic. □ 

Before concluding this section, note that in view of the definition of the function 
Tj^g in (|1.17p . Proposition 11.81 gives the following corollary. 

Corollary 2.2. For a^/p G N and A; G N, there exists a constant C — C'{p, k) such 
that, for a; G R, j G {1, • • • , P} and £ G {1, • • • ,dj}, 

\dfd^r,At,x)\ < e^* (^)(i+P+'=)(i+"o) . 

3. Analysis of the correction terms 
In this section, we will make use of the following norms defined for p G N, 



= sup 



|xre^/(«(x) 



L2 

We associate with this norm the functional space defined by 

= {/ G L2(R'^), \\f\\^.<^}. 
In view of (|1.8p and (11.131) . for all p G N, there exists c{p) such that 

(3.1) ll-^^WIlEf <e^(P)*, yt^O. 

We can obviously take c(0) = by conservation of the L-^-norni, but in general, the 
norm of ip'^ in potentially grows exponentially in time (see [3j). We denote by 
U^{t) the semi-group associated with the operator —^d^ + \k{x) and we observe 
that for p G N, there exists a constant C{p) such that 

(3.2) mmcim^c;ip)e^^p^\'\. 

The following averaging lemma shows an asymptotic orthogonality property. 
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Lemma 3.1. For T > and k ^ j , there exists a constant C such that 



yt£ [0,T], VpeN, 



- f Ul{-s)UI{s)ds 
»t Jo 



^ Ce 



ct 



Proof. We first observe that 

(3.3) ledt {Ul{-t)UI{t)) = [/!(-<) (A,(x) - Xu{x)) [//(t). 

Indeed, if / e L^iK) and f^{t) = UI{-t)U^{t)f. We liave 



iedtrit,x) = - (^--9^ + Afc(x) j r{t) + UU-t) i^--d', + j U^{t)f 

= UI{-t){X,{x)-Xk{x))U^t)f 

because U^{—t) commutes with —^d^ + Xk{x). We use Equation p.3l) to perform 
an integration by parts: 

m{-t)u]it) = ui{-t){x,-Xk)-'umm{-t){x,-Xk)u^{t) 

= ieUl{-t) (A, - Xk)-^ Um dt {m{-t)U^{t)) . 



Therefore, 



1 



Set 
(3.4) 



(A,-A,)-^C/;(.)J^ 
ds (Uli-s) (A, - X,)-' Ulis)) Uli-s)U^{s)ds. 



7j,fe = (Afe-Aj) . 

The behavior as x goes to infinity of these functions is studied in Appendix [C] (see 
Lemma rC.l|) . It is proven there that for all /3 G N, 

Since the propagators U^{t) and Uj{t) map continuously into itself uniformly 
with respect to e, we have 

1 1 



Hi-sh,,kU^is)]l 



els. 



(p+1)tiO+P 



< 



Cip), 



where in all this paragraph, C{p) denotes a generic constant depending only on the 
parameter p G N. Besides, we observe that 



In view of 
1 

ie 

and of 



1 

ie 



= iijAx)ed^+ieil,,{x), 



||f^fc(-s)7j,fe(a;)£<9xC//(s)||^^^(p+3)„o+P+2_j.p^ 
+ ||t^l(-s)7"fc(a^)C^;(s)||£(s(P+2, 



o+p+2 „p^ 
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which comes from p.2[) and Lemma FC-H we get 

\\ds (t/|(-s)7,,fc;7|(s))||^(5..+._5..) < e", 
which concludes the proof. 

We now prove the following proposition. 



□ 



Proposition 3.2. For p e N, there exists C{p) such that for all j ^ 2, and all 

£e{i,...,dj}, 

ll5l,^WIIs?<e^('')*, Vt^O, 
where g| ^ is defined in (|1.16p . 

Proof. We use Duhamel's formula and write 

aliit) = -J u^(t - i'P'isyM) ds. 

it Jq 

Besides, if (pj g{t,x) — Lp^{t,x)rj^i{t,x), then we have, 
[iedt + -dl - Ai(x)j = ledtr,,,^' + r,^,e^'^\^'\^p' + - [dlr,^,{t,x)\ . 



Therefore, we can write 



^U^) = UimiM - I I Ul{t - s)r{s)ds 

whence 

9liit)^-J Unt-s)Ul{s)ds^UQ)--l / W{t-s)Ul{s-T)r{T)dTds 



le 

i / U^^{t^s)Ul{s)ds^lM- 



m(t- s)Ul{s-T)ds 



r^{T)dT. 



Lemma 13.11 yields 



with q = p + 2 + {p + + no). Let us now study r^ . We write = rf + with 

r^(i,a;) = idtrj.iip" + - [dl,rjjXt,x)\ ip" . 
In view of Corollarv 12.21 and of p.ip . we have for all q g N, 



A very rough estimate yields 



< Vie" 



^^^,+(l+,)(„o + l) 



\\{ed.)\' 



where we have used CoroUary[2;i Now with ([113)) and ([3J|) . we conclude 

This completes the proof of Proposition 13.21 □ 
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4. Consistency 

We now prove Theorem 11.91 We go back to Equation (|1.18|) , that we recah: 



.2 

2 

t=0 — 'Oj 



iedte'it^x) + ^—dlem,x) = V{x)e''{t,x) + sNL'^it.x) + eL%t,x), 



where NL^ and are defined in (|1.19p and (ll.20p . respectively. The standard 
L^-estimate yields: 

ii^'wiil^ ^ ii'-oiu^ + f mL'{s)\\L- + w{s)\\L-)ds. 

Jo 

In view of (|1.20p . Proposition 12. II and Proposition 13.21 we have 

\\LHt)U. < Vie^*. 

Besides, we observe 

< (ii^^wiiioo + wmu + e'wfmi^) (iir(t)iu. + s\\fit)u.) . 

In view of (11.131) . we have ||(y5'^(t)||Loo < e~^/^e'-^*. On the other hand, Proposi- 
tion implies, in view of the Gagliardo-Nirenberg inequality 

(4.1) \\f\\L^<e-'/'\\f\\%%d,f\\%\ 
the estimate 

e'\\9%ml^<ee^K 
Therefore, it is natural to perform a bootstrap argument assuming, say 

(4.2) WmL^ ^e-^^^e""'. 

Note that we fixed the value of the constant in factor of the right hand side equal to 
one. We did so because 6*^, as an error term, is expected to be smaller than ip'^ (the 
approximate solution) in the limit e — )■ 0. As long as (|4.2p holds, the L'^-estimate 
implies, in view of (II. 5p 

ll^"(i)llL= /*(Vie^^+e^^||r(s)||i2)ds. 
Jo 

By Gronwall Lemma, we obtain 

(4.3) \\ge^t)h^-^C{e^ + V^)e''''\ 

It remains to check how long the bootstrap assumption (14. 2p holds. For this, we 
use Gagliardo-Nirenberg inequality (14.11) , and we look for a control of the norm of 
O'^it) in Sj. Differentiating the system (ll.lSp with respect to x, we find 

iedtied^e') + -jdl{ed^9') = V{x)ed^e' + eV'{x)e^ + e^d^NL^ + e^d^L', 

We observe that since V is at most quadratic, \V' {x)9^\c;n < (x) |^^|c"- Therefore, 
in order to obtain a closed system of estimates, we consider the equation satisfied 
by xO": muhiply (ITTTS)) by x, 

iedtixO") + —dlixe") = V{x){xe') + e'^dj' + exNL^ + exL". 
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By Proposition [321 'we have 

\\xU{t)\\L^ + \\ed,U{t)\\L^<^ee''K 

Besides, 

\xNL'{t,x)\cM < {\<f>'it,x)\^ + \9'{t,x)\l,^+s^\g'{t,x)\l,^) x 
X {\xe%t,x)\cN +e\xg%t, x)\cn) , 
\ed,NL'{t, x)\cN < + \e'{t,x)\l,^ +e^\g'{t,x)\l,.) \ed,9' {t,x)\cN 

+ s{\<j>'it,x)f + \9'{t,x)\l.+s^\g'{t,x)\l.)\ed^g'it,x)\cN 
+ \ed^(l>'(t,x)\ X X \e'{t,x)\cN 

+ £|0"(i,x)|2 X \d^xHt,x)\c« X \0'{t,x)\cN. 
Arguing as before and using again p.l3p . wc obtain that under ()1.2p we have 

\\ed,e^t)\\L2 + \\xe^t)\\L2 < (e" + V^) e<=^\ 
Gaghardo-Nirenberg inequality then implies 

We infer that (|4.2p holds (at least) as long as 

+ i) e-'" « e-i/^e^*, 

which is ensured provided that t ^ Cloglog (i), for some suitable constant C, since 
K > 1/4. This concludes the bootstrap argument: we infer 

sup {\\d'{t)\\L2 + \\xe'{t)\\L2 + \\ed^9'{t)\\L2)^0. 
ItKciogiog(i) "-^^ 

Theorem 11.91 then follows from the above asymptotics, together with the relation 
= w'^ + eg'^, and Proposition 13.21 

Remark 4.1. In the case where ~ D + eW, as in Remark ll.lOi the proof can 
be adapted, in order to reproduce the argument given in The main point to 
notice is that (local in time) Strichartz estimates are available for the propagator 
associated to —^d^ + D{x), thanks to |6]. Then in the presence of the power e 
in front of W, the potential eW can be considered as a source term in the error 
estimates: the factor e is crucial to avoid a singular power of e due to the presence 
of £ in front of the time derivative in (ll.lSp . The proof in [41 Section 6] for the 
cubic, one-dimensional Schrodinger equation can be reproduced: another bootstrap 
argument can be invoked, which does not involve Gagliardo-Nirenberg inequalities, 
since a useful a priori estimate for the envelope u is available. 

5. Superposition 

As explained in the introduction, the only difficulty in the proof of Theorem ll.lll 
is to treat a nonlinear interaction term. Indeed, we set 

where g'^ is the sum of two correction terms, similar to the one introduced in §1.31 
More precisely, set p{l) = 1, and p{2) = 1 if Ai = A2, p{2) = 2 otherwise. Define 
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f ^ 9i+9h with 

where for k = {1,2}, j ^ p{k) and 1 ^ £ ^ rfj, the function gji,^{t,x) solves the 
scalar Schrodinger equation 



(5.1) i£dt9lkj + Y^^^lkJ- >^j{x)glk,e = ^^rj^kj ; glkj\t=o = 
where 

(5.2) r,-fe^£(i,x) = -i {dtx''it,x)+^p(^k)it)d,x''{t,x) , Xj(i,a:))^, 
The function w'^ (t) then solves 



0, 



w 



|t=0 



0, 



with 



L"^ = ©(Vie' 



E E E 

fc=l,2 isSjXP l^fsjdj 



Here, the ©(-/ee'-^*) holds in S^, from Proposition 13.21 Besides, 

-|<^?|W-|</'2lV2X') 

Adding and subtracting the term ^feYp\'y^ + y'lx^Pl'PiX"'^ + '^'Ix^): we have 

ITVL'^I iV| + iV|, 
where we have the pointwise estimates 

A^l< Vi(|^?lV2l + l'^ll>II), 

N%<^e (l^f p + l^^p + + + £|/|) . 

The semilinear term can be treated exactly in the same manner as in Sectional 

It remains to analyze / ||7Vif (s)|| ds. We observe 
Jo 

x\{s)-xi{sy ^ 



/ ||l'^!(s)lV2(s)||r2ds = 



ds, 



L2 



and we note that the contribution of If/jfPv'l ^^'^ that of I'/slP'/'i play the same 
role. Also, we leave out the other terms which are needed in view of a estimate, 
since they create no trouble. Arguing as in 4^ Lemma 6.1], we obtain: 

Lemma 5.1. Lef T G R, < 7 < 1/2 and 

I%T) ^ {t e [0,T], \x,it)~X2it)\^s-'}. 
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Then, for all k e N, there exists a constant Ck such that 

rT 







with 



Mfe(r) = sup{|| (x)"9f7/j|U»([o^T],L2(R)); je{l,2}, a + P^k). 
In view of this lemma and of Equation (jl.l2p , we obtain 

\\NL\[t)\\^idt < e^^ (^refc(i/2-7) + |r(T)|) , 

and the next lemma yields the conclusion. 
Lemma 5.2. Set 

T = inf \i{x) - X2{x) - {El - E2) , 
and suppose F > 0. Then for < 7 < 1/2, there exists Co, Ci > such that 
\I'{t)\<e'T-'^e^°\ s$ i Cilog f i 



Proof. Consider J^(t) an interval of maximal length included in P{t), and N^(t) 
the number of such intervals. The result comes from the estimate 

|/^(t)| N^{t) X max|J^(i)|, 

with 

(5.3) \r[t)\ < e^e^*r-i and N%t) < te^'T~\ 

provided that e''e^* ^ 1. Let us prove the first property: consider r, cr e J'^{t)- 
There exists t* E [t, a] such that 

|(:^i(r) - x^ir)) - - x^i^))] = \t - <j\ ^t*) - UHl , 

whence 

On the other hand, 
We use 

- = 2 - - Ai(xir)) + A2(a;2(r))) , 

and infer 

Si - ^2 - Ai(:Ei(r)) + A2(a;2(r))| > - ^2 - Ai(xi(r)) + \2{xx{t*)) 

- A2(xi(r))+A2(x2r)) 

^ F - Ce^e^*, 
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where we have used the fact that A2 is at most quadratic. Therefore, if e'^e'-^^ is 
sufficiently small, 



El - E2 -\i{xi{t*)) +M{x2{t*)) 



r 

^ -. 
2 



We infer 

<e''e'^*r-i, 

provided e^e*"* 1. 

Let us now consider N^{t). We use that as t is large, N^{t) is comparable to 
the number of distinct intervals of maximal size where \xi{t) — X2{t)\ ^ e^. More 
precisely, N'^^t) is smaller than t divided by the minimal size of these intervals. 
Therefore, we consider one interval ]t, cr[ of this type and we look for lower bound 
oi a — T. We have 

|xi(r) - X2(r)| = |xi((t) - a;2(cr)| = e''', and Vt G [r, cr] , \xi{t) - X2{t)\ ^ e"' . 

Besides, inside ]t, a[, xi{t) — X2{t) has a constant sign that we can suppose to be + 
(one argues similarly if it is — ). Under this assumption, we have 

- > and Ua) - < 0. 

Using the exponential control of X'j{xj{t)) for j e {1, 2}, we obtain 

(5.4) (6(r) - 6(r)) - (6 (a) - U<^)) < e«(a - r). 

We write 



|l6(r)p-|6(r)pl 



(5.5) a(r)-6(T) = |ei(r)-6(r)| > 



ia(r)| + |6(r)| 
>e-^*||a(T)p-|6(r)n 

and 

(5.6) - 6(^) +6('t) = Mr) - 6(t)| > e-" ||6('t)P - 16(^)1'! • 
As before, we prove 

|l^i(r)p - I6(r)n + |ICi(^)P - I6(^)n > r, 

provided that e'''e'"* <C 1. Therefore, plugging the latter equation, (|5.5p and (15.6 
into (15.41). we obtain 



CT - r > e-2Ctr thus N'{t) < te'^'T'^ < e^*r~^ 
which completes the proof of Theorem 11.111 □ 



Remark 5.3. The proof shows that if the approximation of Theorem 11.91 is proven 
to be valid on some time interval [0, Clog(l/e)], then Theorem 11.111 will also be 
valid on a time interval of the same form. 
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Appendix A. Global existence of the exact solution 

The proof of Lemma [1.31 follows classical arguments; see [15] (or [5]) for more 
details. We suppose e = 1 without loss of generality. We use the decomposi- 
tion V{x) = D{x) + W{x) of Assumption 11.21 and we denote by U{t) the unitary 
propagator of —^9^ + D{x). Let Xt be the set 

Xt = {V' e C(/t,S}), V,a;V',VV' e L«(/T,i''(R,C^))} , It ^]s-T,s + T[ 

for s G R and T G R to be fixed later. The proof consists in a fixed point argument 
for the function 

where for s G R, the function $s('(/') is defined by 



^s{-il^){t) = U{t-s)^{s)-iK / U{t-T){\i:\},«^){T)dT-l / U{t~T){W^{T))dT. 



By [6], local in time Strichartz estimates are available for U. Strichartz estimates 
and Holder inequality imply that there exists a constant C > such that 

l|4's(V')l|LS(/T,-L*)nL-(/T,i^) «S C\\iI}{s)\\l2 + C||?A||^s/3(/^,i4)||V'||LS(/r,L'') 

+ C\\WnLHi^,L-y 

Using the boundedness of the coefficients of W and Holder inequality in time, we 
obtain 



We can then infer that is a contraction on a ball of Xt for some T which depends 
only on |lV'(s)||i2. Then, the conservation of ||-0(t)||2,2 yields the lemma. 

Appendix B. Some formulas involving the projectors 

In this section, we list and prove some formulas which will be used in the course 
of the computations in the next appendix. We consider here the more general case 
X € R'^, with d ^ 1. Fix once and for all in this paragraph j G {1, . ■ ■ ,P} and 
i G {1, . . . ,d}. First, recall that we have seen in ^that since itj = Hj, 

(B.l) H,(9,n,)H, =0. 

Differentiating the relation H| = Uj, we find: Vj G {1, . . . , F}, G {1, . . . , d}, 
(B.2) deU, = (deU, + H, (9,H, ) . 

We now show: Vj G {1, . . . , P}, G {1, . . . , d}, 



= (nfe(9,n,)H,+n,(a,n,)Hfe), 

l<k<P 



(B.3) 



where the last equality stems from (|B.1|) . To prove (jB.Sp . simply write 

a,H, = ^Hfe(a,nj)H„ 

k.m 

where we have used J^k ^fe = Id- Then, observe that HfcHj = Sj^Hj yields 
Ukidell^) + idink)Iij whence UkideUj) = -{diUk)Uj. 
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The fact that Hj = for all m ^ j gives (jBTSj) . 

The last formulas we wish to establish involve the spectral gap. Since we have a 
basis of eigenfunctions, we have 

VU,=Il,V = X,Il,. 

Differentiating with respect to xe, we infer 

idenj)V + IljdeV ^ XjOiR^ + {diXj)nj. 

For k S {1, . . . , P}, multiply this relation by Hk on the right, and use the property 

Afc(a,n,)nfc + n, {deV - dAj) = Aj(a,n,)nfe, 

hence 

(B.4) (A, - Xk){deIlj)Uk = {diV - diXj)Uk. 

Similarly, we have 

(B.5) (A, - Xk)Tlk{den,) = Uk (diV - 9, A,) n,. 

Appendix C. About the growth of the eigenvectors at infinity 

This section is devoted to the proof of estimates at infinity for the eigenprojectors 
associated with a potential V satisfying Assumption 11.21 We will use a lemma on 
the derivatives of the inverse of the gap between two different eigenvalues. For 
j, fc G {1, . . . , P}, j k, we recall that we have set (see (13. 4p ') 

Vx g R, ^j^kix) ^ {Xj{x) - Xk{x))~^. 

Since the results are not specific to the space dimension one, we prove them for 
potentials depending on a; e R*^, d ^ 1. 



Lemma C.l. Assume (jl.2p is satisfied with rip G N and that the functions V 
and Xj (j £ {!, - ■ ■ ,P}) are at most quadratic. Then, for /3 G N'' and for j, k G 
{1,...,P} with j ^ k, 

(C.l) \d!^lj,k{^)\ ^ (^)«o + l^l(l+no) _ 

Proof. For /3 = 1^, we immediately obtain 

di (Xjix) - Xk{x)) 



\9ajA^)\ 



iX,{x) - Xkix)f 



from (|1.2p . and the fact that Xj and A^ are at most quadratic. 

Set Aj^k = Xj — Xk- it is at most quadratic. Besides, for (3 G N'', we have 



c»l + --- + Qp=/3 



for some real-numbers aaj^^,,,^ap- The result then follows by observing that 

ATi-^'a- A,,, . . . a:^A,,,| < (.)-(^+^) {xr , 

from (|1.2p . and the property |i9"Aj^fc| < (a;)*'^"'"'''*, which follows from the fact 
that Aj^k is at most quadratic, in the sense of Definition ll.il □ 
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We now consider the eigenprojectors IIj associated with the eigenvalues \j of 
the matrix V. Because of the gap condition, these functions are smooth in R'^. We 
prove the following 

Lemma C.2. Let Hj be an eigenprojector of V for j € {!,..., P}, we have for 

(C.2) |5fn|c«,« < (x) 1^1(1+""), 

where the norm \ ■ |c"-" denotes the matricial norm. 

Proof. The case \f3\ = is immediate since IIj is a projector. In view of (jB.Sp . 
relations ((B31) and (|R4)) imply (1021) for |^| = 1. 

We now argue by induction. We suppose that (|C.2p holds for any 7 e N** 
with I7I = K for some K ^ N and we consider /3 with |/3| = K + 1 and /3e ^ 0. 
Differentiation of order 13 — li of (|B.2[) and multiplication on both sides by IIj yields 



\0<|a|<l/J| / 

where all along this proof, Oa will denote real numbers whose exact value is unim- 
portant. We obtain 

(C.3) |n,(9fn,-)n,lc„.„ <(x)i'^i(^+"°). 

Then, for all k ^ j, we estimate {d^IIj)Ilk. To do so, we differentiate (|B.4|) and get 

0<|a|<l/3 

+ E ^",,...,..9^^ ((A, - dTn,d:^d, (v - A,) 9^ Bfc. 

QlH |-Q4=/3-lf 

In the first sum, the induction assumption yields 

(C.4) ||a^n,c)f-"nfe||^„.„ < (^)l"l(l+"o) + l,9-a|(l+no) ^ ^^^|/3|(l+no) 

Besides, for each term in the second sum, we write 

11^:^ ((A, - Afe)-i) arn,9,"3a^^ „ A,) 9rnfc||c«,« 

< /^^'■^no + \ai\(l+no) ^^\j(l-|a3|)+ (1+no) (|a2 1 + |a4 1 ) 

where ?'+ = max(r, 0) and where we have used the fact that V and Xj are at most 
quadratic, together with the induction assumption and Lemma IC. 11 We have the 
two alternatives: 

• If q;3 = 0, then 

no + |ai|(l + no) + (1 - |a3|)+ + (1 + no)(|a2| + |a4|) 

= no + |ai|(l + no) + 1 + (1 + no)(|a2| + |a4|) 
= (l + no)(l + |ai| + |a2| + |a3|) = (l + no)|/3|, 
since ai + a2 + a4 — (3 ~ l£. 
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• If as 7^ 0, then 

no + |ai|(l + no) + (1 - |a3|)+ + (1 + no)(|a2| + |a4|) 

= no + (1 + no)(|Q;i| + \a2\ + |a4|) 
^ (l + no)(l + |ai| + |a2| + |a4|) 
^ (l + no)|/3|. 

We deduce 

(C.5) ||(9fn,)n,.||^„,„ < , Vfc ^ J. 

Similarly, 

(C.6) ||nfe(5f n,)||^„,„ < , Vfc ^ J. 

In view of (|C.3p . we infer 

(C.7) ||n,(afn,)nfe||^„.„ + ||n,(afn,)n,||^„,„ < , yj^k. 

Applying the operator d^~^^ to (|B.3|) . the induction assumption and equations 
((U31) . ([CJ]) and ((CJ)) yield ((a2|) . which concludes the induction. □ 
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